Self-Similar Solutions to a Density-Dependent Reaction-Diffusion Model 
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In this paper, we investigated a density-dependent reaction-diffusion equation, ut = (u m ) xx + 
u — u m . This equation is known as the extension of the Fisher or Kolmogoroff-Petrovsky-Piscounoff 
equation which is widely used in the population dynamics, combustion theory and plasma physics. 
By employing the suitable transformation, this equation was mapped to the anomalous diffusion 
equation where the nonlinear reaction term was eliminated. Due to its simpler form, some exact 
self-similar solutions with the compact support have been obtained. The solutions, evolving from 
an initial state, converge to the usual traveling wave at a certain transition time. Hence, it is quite 
clear the connection between the self-similar solution and the traveling wave solution from these 
results. Moreover, the solutions were found in the manner that either propagates to the right or 
propagates to the left. Furthermore, the two solutions form a symmetric solution, expanding in 
both directions. The application on the spatiotemporal pattern formation in biological population 
has been mainly focused. 

PACS numbers: 82.40.Ck, 87.23.Cc, 02.30.Jr, 87.18. Hf 



I. INTRODUCTION 

Many transport phenomena arising in natural sci- 
ence have been successfully described by the nonlinear 
reaction-diffusion (NRD) model. The applications of the 
NRD model have had wide variety including population 
dynamics, transport in porous medium, combustion the- 
ory and plasma physics The NRD equation is for- 
mulated by 



ut = [D(u)u x ] x + R(u), 



(1) 



where u :— u(x,t) ^ 0, depending on the system, is den- 
sity or temperature at position x G (— oo, oo) and time 
t £ [0, oo). D(u) and R(u) are the diffusion coefficient 
and the reaction term respectively. The simplest form of 
NRD Eq. |T]) is the well known Fisher equation Q , as well 
as the Kolmogoroff-Petrovsky-Piscounoff (KPP) equa- 
tion, which was originally derived for the propagation of 
a gene in a population. In this classic model, the reaction 
term is represented by the self-limiting Pearl- Verhulst lo- 
gistic law R(u) — au(~l — u/a) and the diffusion is rep- 
resented by the Brownian process, reflected through the 
term D(u) = k, where k > is the diffusion constant, 
a > is rate constant and a = lim u(x, t) is equilibrium 

t— ► oo 

density. Since the Brownian motion seems to be not real- 
istic movement for the active individual such as biologi- 
cal population, the density-dependent diffusion has been 
included which represented by D(u) = k(u/o-) p where 
p > 04H| ■ It reflects the crowd avoiding movement of 
the individual. Later, it has been also found the more 
general form for the logistic law R(u) — au[l — (u/a) p ] 
as well [llj]. By substituting D(u) and R(u) in the 



later case into Eq. ([T]) with the transformations t* = at, 
x* = (ma/ k)^ x and u* = u/a, we obtain the dimension- 
less equation 



u t = (u m ) xx + 



(2) 
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where m = p + 1 > 1 and the asterisk is dropped for con- 
venience. The applications of this equation have been 
found in various systems (see [U-IH for instance). Due 
to the diffusion coefficient is depended on the density; 
this equation is so called the density-dependent reaction- 
diffusion (DDRD) equation. Sometime, this model is 
called the generalized or extended Fisher equation Q. 
Since the existence of the traveling wave solution in the 
NRD equation has been found in both the original Fisher 
equation [(| and the DDRD equation [lol [TlT | , the analy- 
sis on these equations has been attracted much attention 
(see [1, [|[ for review) . 

Recently, the NRD model has been employed to de- 
scribe the spatiotemporal pattern formation in bacterial 
colonies which is the ele gan t self-organization observed 
in the biological system |12h14| . For each species, the 
bacterial colonies generate various patterns depending on 
surrounding environments, mainly, including the nutrient 
concentration and the softness of the medium. This coop- 
erative behavior in adaptation to survival in the environ- 
ment reflects bacterial communication capabilities and 
social intelligence [TH, [n| . The understanding of the un- 
derlying mechanism is not only important to the biotech- 
nology but also the basic science of the living organisms. 
Kawasaki et al. [HJ pointed out that although the bac- 
terial colony evolves in two dimensions, each tip grows in 
one dimension except for occasional branching. There- 
fore, the velocity of the tip elongation obtained form the 
one-dimensional system could be compared well with one 
obtained from the two-dimensional system. They also 
proposed the simplified version of the model for pattern 
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formation in bacterial colony which was consistent with 
Eq. ©0,0. 

In this paper, we investigate the generalized or ex- 
tended Fisher equation (Eq. ([2])) in one-dimensional 
space [H [ll|. Although the asymptotic behavior of its 
solution has well understood as the traveling wave solu- 
tions from the past studies, the exact or event explicit 
solution of Eq. ([2]) has been unknown. That hinders the 
understanding how the systems develop to the traveling 
wave from the initial state. In different aspect, the self- 
similar solution has been substituted directly to Eq. © 
16] . It turned out that, after evaluating, the solution was 
difficult to get the closed-form. The analysis on the re- 
duced form of Eq. © has been also carried out. In 1977, 
Gurtin and MacCamy Q have shown that the NRD, with 
Malthusian law as the reaction term, can be reduced to 
the simpler form by the simple transformation, which the 
reaction term was eliminated. Thus, equation transforms 
to the purely nonlinear diffusion equation. In 1983, New- 
man [ljj has transformed Eq. © to the reduced form un- 
til the integral invariant and Lyapunov functional have 
been admissible. The convergence of the solutions toward 
to the traveling wave has been obtained as the asymp- 
totic behavior. In 2002, Rosenau [l?} has employed a 
transformation which mapped Eq. © to the purely non- 
linear diffusion where the the nonlinear reaction term 
R(u) was eliminated. According to the inhomogeneous 
nature of the mapped equation, the solution has been 
obtained only at asymptotic form and has not been in 
the closed-form. In 2004, Harris 18| has also applied 
the similar method of Rosenau [l7 to analyze Eq. ©. 
Indeed, the explicit solution, which evolves from the ini- 
tial state to the traveling wave, has been obtained for 
to = 2. However, the solution for the case to > 2 has 
been unknown. 

In our analysis, we also follow the similar method of 
[i~7l Qjj]. But, the Eq. © is differently mapped to the 
anomalous diffusion equation which the exact self-similar 
solution is obtained for any m > 2 (l^ - [2lj |. It is im- 
portant to point out that the connection between the 
self-similar solution to the traveling wave solution, as de- 
scribed in [l^], has seen clearly from this particular so- 
lution. The outlines of this paper are as followings. The 
details of analysis for the solution of Eq. © are described 
in Sec. |Hj The discussions on the obtained results are 
given in Sec. IIII1 Finally, we conclude our study in Sec. 

Wl 



II. MAIN RESULTS 



We now perform the analysis on Eq. ©. Noticing that, 
Eq. © can be factorized as (^-l)it = (^ + 1)( J^-l)u m 
and then it can be written in the equivalent equation 



dt 



. d x . x d _ 
dx dx 



(3) 



It suggests that Eq. © can be reduced to the simpler 
form by seeking the solution in the composition of the 
exponential factors 



u(x, t) = e t e™$(x, t), 



(4) 



where $(a;,t) ^ is the transformed density. By substi- 
tuting Eq. © into Eq. ([3]), we obtain 



e -(m-l)t $ 



- X [e 2x {^ m ) x ] x . 



(5) 



To eliminate the exponential terms in Eq. (J5]) , we define 
(to — l)e^ m ~ 1 )*di = dr and then the temporal transform- 
function is given by 



r(t) 



a (m-l)t 



1. 



(6) 



which satisfies the initial condition r(0) = 0. Similarly, 
we define ^- L ^-e^~ x dx = dd> and then we evaluate the 
spatial transform-function 



4>{x) 



(7) 



With transformations Eq. ((6]) and Eq. ©, the intrinsic 
space and time are eliminated. Consequently, Eq. © is 
reduced to 



r 771 — 1 



m+l > 



[0^r($™),] . 



(8) 



Eq. © is similar to the anomalous diffusion equation 
$ r = [fc^ 1 where k = (™±±) 2 and I = 

> 2 ' which the exact solution is known [l9l - [2l| . To 
this end, it can be interpreted that the nonlinear reaction- 
diffusion Eq. © is mapped to the anomalous diffusion 
Eq. © by employing the transformations Eq. ©, Eq. © 
and Eq. ©. 

Similar to (I9I - I21I ] , the solution of Eq. © is assumed 
to be the normalized scaled function of the type 



1 



T{tY \T(t)J T{t) 



-F 



(9) 



where uj((f>, r) :— 0/T(r). By substituting Eq. © into 
Eq. ©, we obtain 



rp2 y u ' m— 1 V m 



n a [oj^r(F m ) u 



rp2+r. 



(10) 



Eq. (|10l) can be separated into two equations; The first 
and the second equations are given by 



T m[m—X) 



(11) 



and 



X^^^F^^^Lo^+riF™)^, (12) 



respectively, where A is an arbitrary positive constant. 
According to Eq. (fTTj) , we obtain the solution 



T(r) = X(T + a)^-i) 



(13) 
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where a is an arbitrary positive constant. By integrating 
both sides of Eq. (TT2")) . we obtain 



Tn(m-l) . i 3in + l 



A 



E±± w W(F m ) |l 

1) 



(14) 



which is equivalent to — ( 
a; — > — oo, we have w(0,t) = <p(x)/T(r 



0. From 



Eq. (|T4| . the boundary condition is shown F w (0) = 0. 
Again, Eq. (IT4"|) can be solved for the solution 



F(w) = A b+ (Acj)^^tt 



(15) 



where 6 is a constant which is determined later. By sub- 
stituting Eq. flT3J) and Eq. dT5j) into Eq. ©, we obtain 
the transformed density 



3(6 t) 



1 



(r + a)- 



(r + a)" 



(16) 



We substitute Eq. (fT6|) into the main solution Eq. Q , by 
doing some algebra, then we obtain 



u(x, t) 



(ef* - 1 + a)p 



1 



(17) 

which is the implicit solution of Eq. . 

Our selection of the undetermined constants a and 
b are dependent on the initial- and boundary con- 
ditions. At t = 0, we have the initial profile as 



1 



■b(eP x /a) 1/{p+1) 



If 



u {x) := u{x,0) = (l/a 1 ^) 

the solution has the (semi) compact support 
the density must be vanish at an initial front posi- 
tion xo, i.e. uq(xq) = 0. We then evaluate the con- 
stant b = -(ae-P X0 ) 1 ^ ( P +1 \ The initial density now 
becomes M (x) = (1/a) 1 ^ [l - e P(*-*o)/(p+i)] As 
x <g; xo, density must reach the constant value p = 
lim uo(x) — (l/a) 1 / p . That is, we have a = p~ p . 

By substituting the constant a and b into Eq. (|TT)) . we 
obtain the exact solution 



u(x, i) 



1 - 



pe 



[pP(eP* - 1) + 1] 

„p(a;-a; ) 



_pP( e P* - 1) + 1_ 

which the initial density profile has the form 

i 

*r(z-sa)l p 



u (x) = p 



1 



(18) 



(19) 



As x <C xq , the density converges to the equilibrium state 
lim u(x, t) = 1. 

t— »oc 

We have also found another solution for Eq. ^ . In the 
factorized form, we rearrange Eq. ^ as — l)lt = (^j — 



*0 -t 



1)(^: + l)u m . Then, it can be written in the equivalent 
equation 

d d 

6 Tx e (e & e " } ' (20) 

By omitting the details, similarly, by defining another 
solution u_(x,t) = e~™e'$(x,i) for Eq. (0), we ob- 
tained the same anomalous diffusion equation Eq. ([8]) and 
the same temporal transform-function Eq. (|6]) with the 
different spatial transform-function 0_(x) = e — 
4>(—x). Later, it is found the relation between the two 
solutions actually that 



- (x, t) = u{— x, t). 



(21) 



It is obvious that the two solutions are mutual sym- 
metry. Once, it was sugges ted that we can construct the 
16l. \vf\ by the linear combination 



symmetric solution 1^, 
of 



w(x,t) :=u{x,t)+u-{x,t), (22) 

To see how Eq. (|2"2"|) applicable is, let us consider the 
combination of Eq. @ of the two solutions (it + U-)t — 
(u m +u™) xx + (u+u-) — (u m +u™). By using the binomial 



expansion u 



for any integer m and u 

Eoo jm 
fe=l \ k 

the approximation 



= («+tt_r-E^ (?)» 



-k„,k 



= (u + u_) m - 



u m k u k _ — tt™] for any fraction to, we obtain 



u m + u m w (u + u-) m + 0{H), (23) 

where H(x,t) = (™)u fe (x, t)it™~ fc (x, i) for integer 

to and H(x,t) = J2T=i h^'Hx,^^'^ - u™(x,t) 
for fraction to. Using Eq. (|23[) . then we obtain 



I0 t = K) XI + 



10 



0{H xx -H). (24) 



Therefore, the linear combination w(x,t) in Eq. (|22|) can 
be the solution of Eq. © with the correction 0(H xx — H). 
For any p > 0, using the binomial expansion Eq. (|23l) and 
modifying the initial condition, we have the symmetric 
solution 

np t 

w(x, t) 



pe 1 



1 - 



[pp( e pt - 1) + i]p 

^^-(X-Xp) 



[pP( e P* _ 1) 
= 4(6*+* 



(25) 



f l] p+1 J 

where Ch(x + k) := ^(e x+k + e~ x+k ) = e fe coshx is the 
modified hyperbolic function. It is found that the sym- 
metric initial density uiq(x) := w(x,0) is given by 



w (x) 



1 - Ch- 



P 



P 



1 



(26) 



Here, we found the special forms of the symmetric initial 
density u>o(x). If xo is sufficiently small io < 1, the 
initial density forms the delta- like profile wq(x) fts p5(x). 
If xo and p are sufficiently large, xo 3> 1 and p > 1, 
the initial density forms the rectangular profile Wq(x) k, 

P n{x). 
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FIG. 1. The spatiotemporal evolution of the density profile in 
the case p = 4 with the initial conditions p — 0.2 and xa — 1. 
The solid lines represent u(x, t) and the dashed lines represent 
u- (x, t) — u(—x, t). 




FIG. 2. The spreading velocity v(t) of u(x,t) corresponding 
to Fig. HJ. The initial speed is v(0) = p 4 = 0.0016 and the 
transition time is t' w — lnp = 1.61. 



III. DISCUSSIONS 

The evolution of solution u(a;, t) in Eq. (fl8|) and its cor- 
responding U-(x,t) = u(—X,t) are illustrated in Fig. ([1}. 
It has seen that u(x, t) is spreading to the right whereas 
u~(x,t) is spreading to the left, with velocity iv(t) re- 
spectively. The solutions grow to the saturated density 
profile with the maximum value 1 by initiating p < 1 . In 
other hand, the solutions decay to the saturated density 
profile with the maximum value 1 by initiating p > 1 
(has not been illustrated in the figure here). 



A. Spreading velocity and transition time 



As the solution Eq. (|18|) has the compact support, the 
density falls to zero u(r, t) = after the front position 
r(t). Consequently, we obtain the front position for the 
right-spreading solution only: 



r(t) = x + 



\n(pP(e pt - 1) + 1) 



(27) 



Then, the spreading velocity is calculated as 
dr(t) 



«(*) = 



p p e pt 



dt 



pP(eP t -!) + !' 



(28) 



The spreading velocity versus time is illustrated in 
Fig. ©. If the initial density constant is not equal to 
the equilibrium point p ^ 1, it is seen that the dynamics 
of the systems can be separated at least in two regimes. 
At the early stage, the systems spread rapidly from an 
initial speed v (0) = p p . After the transition time if , the 
systems spread with the constant the spreading speed 
c = lim v(t) = 1 at the late stage. Therefore, the transi- 

t— >oo 

tion time t' can be estimated from Eq. (|2"5)l by using the 
fact p p {e pt ' - 1) w 1. Then, we obtain 



if 



ln(l 



1 



(29) 



which is exactly depended on the initial density constant 
p. If the initial density constant is much greater than the 
equilibrium density p> 1, the transition time is expected 
t' ps 0. Otherwise, if the initial density constant is much 
lesser than the equilibrium density p 1, transition time 
is expected if ps — In p. 



B. Traveling wave solutions 

At the large time scale t 3> t', the solution Eq. ((T 
becomes 



u{x — ct) 



(30) 



which is known as the traveling wave solution where c = 1 
is the front speed. When converting to the physical di- 
mension, the front speed obtained here c = y/ kol /(p + 1 ) 
is comparable to the minimum value that the traveling 
wave emerged [3, 0, E3, El • If the initial density constant 
is equal to the equilibrium density p = 1, the traveling 
wave emerges from t = without the transition time. 
Similarly, at the large time scale t ^> t' , Eq. (l25t devel- 
ops to the symmetric traveling wave solution 



w(x — ct) 



Cb -^+l( x -t-x ) 



p> 



(31) 



which is the compact expanding wave with the expanding 
speed c = 1. We now have the traveling wave solutions 
either spreading to the right or to the left and the sym- 
metric traveling wave solution spreading both direction, 
as shown by [23j . 

We point out that, at large time scale t^>if, the so- 
lution Eq. (|16[) forms the self-similar structure $(</>, t) ps 
■\F(At), where B — T +1 -, \ , which can be classified as 

T p \ T p/j r m(m — 1) ' 

the intermediate asymptotics of the second kind as de- 
scribed in 0, [2^1 • The traveling wave solution Eq. (|3"0)l 
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u{x — t — xq) is connected to the self-similar solution by 
the transformations u — > e x / m e t $>, r — > e (m-i)t^ anc [ 
6 — > e ( m + 1 )( x ~ x o)/"\ 

C. Biological population dynamics 

In the term of biological system, the solution Eq. (|25|) 
captures the qualitative features of the spatiotemporal 
pattern formation bacterial colonies 14j |. The bacte- 
ria are inoculated from an initial density profile then the 
bacteria reproduce until reach the saturated number due 
to the limitation of nutrient while the bacteria expand 
the colony to localize the higher nutrient region. Surpris- 
ingly, this process has reflected the self-similar structure. 

In the opposite of the logistic law, a < 0, Eq. ([2]) 
becomes Ut = (u m ) xx + u m — u, in the dimensionless 
form. This equation is exactly the porous medium equa- 
tion with source and sink term. We have found that by 
the transformations t —t and x — > ix, absolutely, the 
symmetric solution Eq. (p25)l would be applicable. Since 
the traveling wave pattern has not emerged in this sys- 



tem, the further analysis on this equation is not scope of 
this paper. 

IV. CONCLUSIONS 

We studied the spatiotemporal pattern formation in 
a density-dependent reaction-diffusion equation which is 
known as the extension of the Fisher or Kolmogoroff- 
Petrovsky-Piscounoff equation. The exact self-similar so- 
lutions with the compact support were found. The so- 
lutions converge to the well-known traveling wave solu- 
tions at t — > oo. Therefore, the connection between the 
self-similar solutions and the traveling wave solutions has 
been shown clearly from the analytical solutions. 
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